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We investigated domain growth in switching processes between the low-spin and high-spin phases 
in thermally induced hysteresis loops of spin-crossover (SC) solids. Elastic interactions among 
the molecules induce effective long-range interactions, and thus the boundary condition plays a 
significant role in the dynamics. In contrast to SC systems with periodic boundary conditions, where 
uniform configurations are maintained during the switching process, we found that domain structures 
appear with open boundary conditions. Unlike Ising-like models with short-range interactions, 
domains always grow from the corners of the system. The present clustering mechanism provides 
an insight into the switching dynamics of SC solids, in particular, in nano-scale systems. 

PACS numbers: 75.30.Wx 75.50.Xx 75.60.-d 64.60.-i 




Spin-crossover (SC) compounds have been studied in- 
tensively because of their peculiar physical properties due 
to competition between the low energy of the low-spin 
(LS) state and the high entropy of the high-spin (HS) 
state [EH SC transitions are induced by changes in 
temperature, pressure, etc. The LS state can be excited 
by photo-irradiation to a long-lived HS state at low tem- 
peratures, which is called LIESST (light induced excited 
spin state trapping) and reverse LIESST (HS to LS) 
can also be obtained at a different wavelength 
controllable and functional properties [H, 0, H, 
would bring potential applicability to novel optical de- 
vices, e.g., optical data storage and optical sensors. 

The LS and HS states couple through a vibronic mech- 
anism and the size of the SC molecule changes with 
the spin state. The distortion caused by the change of 
molecular size induces a kind of elastic interaction among 
the spin states of molecules. The importance of the 
elastic interaction has been reported for SC transitions 
0, El El EI 0, 0, 0, 0, El, 0, EH . The mechanism 
of the phase transition induced by the elastic interaction 
has been studied and various new aspects have been re- 
vealed. For example, using periodic boundary conditions, 
it has been shown that effective long-range interactions 
suppress domain growth, and uniform configurations are 
maintained even near the critical temperature [16[. In 
the process of switching the configuration uniformity is 
also maintained, which is considered an intrinsic property 
of systems with the effective long-range interactions. In 
these systems, we do not expect the critical opalescence 
due to growth of large clusters. 

Nowadays SC compounds are a focus of nano- science 
and technology [H, 0, 0, 0, 0]- On the nano-scale, 



such as powder or thin film samples, the boundary effect 
is important. In particular, in systems with long-range 
interactions the concept of the thermodynamic limit may 
not be well defined and the effect of the boundary must 
be considered carefully. 

In this Letter we investigate how a SC system with 
effective long-range interactions switches between the 
bistable states, using open boundary conditions. We an- 
alyze characteristic features of the heating and cooling 
processes in thermal hysteresis loops with open bound- 
ary conditions (OBC), and compare to periodic boundary 
conditions (PBC). 

We adopt a simple SC model for the square lattice, 
which represents general characteristics. In the model, 
both intramolecular and intermolecular interactions are 
taken into account 201, 
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Here, Xi and Pi represent the coordinate and its con- 
jugate momentum of the center of mass for the ith 
molecule. Conjugate variables and Pi are defined 
for the totally symmetric mode for the ith molecule, 
which is the most important intramolecular motion [27]. 
We also define the variable x as x = r — tls, where 
^ls(= 9) is the ideal radius of the LS molecule. That 
of the HS molecule is ths = ^ls + 1- The intramolec- 
ular potential energy ^ mtra (xi) is shown by the solid 
curve in Fig. [I] (a). We adopt the intermolecular poten- 
tial Vij lter ( Xi, Xj,ri, rj) between the nearest neighbors 




FIG. 1: (color online) (a) Intramolecular potential energy 
V(x) shown by the solid (blue) curve. A realistic value 
^| = 2 is adopted [27[. The dotted curves are LS and HS 
potential energies without quantum mixing. A LS molecule 
(blue small circle) and HS molecule (red large circle) are inset, 
(b) Thermal hysteresis loop for D = 20 and L = 100 with the 
open boundary. Open circles denote HS fraction. 



and next-nearest neighbors, where D is the strength of 
the intermolecular interaction [281 ]. 

We study the present model by a molecular dynamics 
method, in which we introduce a mechanism to control 
the large entropy difference between the HS and LS states 
[20[. When the intermolecular interaction is stronger, the 
system exhibits a thermal hysteresis. Using PBC, uni- 
form configurations are maintained during the transition 
between the HS and LS phases. Using OBC, however, a 
macroscopic inhomogeneity is produced. 

Here we focus on the dynamics of a relatively large 
hysteresis loop (D =20). In Fig. [T] (b) the temperature 
dependence of the HS fraction [20[ is given for a sys- 
tem of N = L 2 = 100 x 100. The system was heated 
from T = 0.1 to 3.0 in steps of 0.1, and then cooled to 
the initial temperature T = 0.1. At each temperature, 
first 400,00 MD steps were discarded as transient time, 
and then 200,00 MD steps were used to measure physical 
quantities with the MD time step At = 0.01. The tran- 
sition from the LS to HS state (LS — > HS) occurs around 
T = 1.0 in the heating process and from the HS to LS 
state (HS — > LS) around T = 0.4 in the cooling process. 

We clearly observed domain growth from corners dur- 
ing the transition from the HS to LS phase (Fig. [2] (a)- 
(d)). LS domains grow and finally combine together to 
extend to the whole system. Domain growth occurs in 
the diagonal directions. We studied the system-size de- 
pendence of the process and found that clusters always 
grow from the corners, not from the sides (edges) or the 
inner part (bulk). The qualitatively same characteris- 
tics of clustering of LS domains were observed in relax- 
ation processes from the metastable HS phase at a low 
temperature (T < 0.4), which is realized by LIESST in 
experiments. 

The feature of clustering presents a distinct contrast 
to the cases of Ising-like models with short-range inter- 
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FIG. 2: (color online) Snapshots of configurations in the 
cooling process from the HS to LS phase. Red (blue) circles 
denote HS (LS) molecules, (a) T = 0.4, t = 33108, (b) T = 
0.4, t = 33506, (c) T = 0.3, t = 33692, (d) T = 0.3, t = 33804. 

actions, where the nucleation occurs from the inner part 
or the sides in large systems |29j . In order to clarify this 
difference, we study energy dependence on the cluster 
pattern, and we find that the growth from sides is not 
acceptable in the elastic model. 

As the initial state we set a round LS domain (closed 
quadrant) at a corner in the complete HS phase. Then, 
we move all molecules slowly by reducing the total poten- 
tial energy of the system. We define N^ s as the number 
of LS molecules in the LS domain in the initial state, and 
TV^s as that in the stationary state. R\^ s is the number 
of LS molecules in the horizontal (X) (or vertical (Y)) 
direction for the domain with A^g. 

Figure[3](a) illustrates the configuration in the station- 
ary state when we set = 7 and N*[ s = 39. In this 
case, the quadrant shape was maintained and the number 
of LS molecules in the domain did not change during the 
simulation, i.e., N^ s = N^ s , although a small shift of the 
position and radius for each LS molecule was observed, 
which means that the LS domain in the stationary state 
is at least locally stable. 

We investigate AE, defined as the difference of energies 
between the stationary state and the complete HS phase. 
The dependence of AE on is plotted by circles in 
Fig.[U We find that AE becomes lower as N^ s increases, 
which indicates that the larger LS domain at the corner 
is energetically favorable. We checked that it holds true 
for larger system size (L). Thus, the domain will grow if 
some small noise (thermal fluctuation) assists the system 
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FIG. 3: (color online) (a) Stationary configuration of a round 
LS domain at a corner in HS phase, (b) Initial configura- 
tion of a round LS domain at a side (c) A configuration in 
the intermediate state. Green color denotes unstable state of 
molecules, (d) Configuration in the stationary state. Unsta- 
ble LS molecules changed to HS molecules and the initially 
round LS domain does not exist any more. 
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FIG. 4: (color online) AE vs. iV^s f° r LS domains at the 
corner and the side. For LS domains at the corner, iV^s is the 
same as A^g, unlike LS domains at the side. The number in 
a parenthesis is the value of N^ s when the initial LS domain 
shape is not maintained and N^ s varies from iV£ s . 



to relax. 

Next, we study the stability of LS domains (closed 
semicircle) at a side. Figure [3] (b) shows the initial 
state of the configuration of a LS domain (i?Ls = 7, 
TVls = 78) at a side. In this case, the initial config- 
uration is found to be unstable due to high distortion. 
Several LS molecules change back to HS molecules and 
the number of LS molecules changes (iV^g > N^ s ). 

Figure [3] (c) is a configuration in the intermediate state 
and the final stable configuration is shown in Fig. [3] (d). 
The molecules colored green (in gray in black-and-white 
print) in Fig. [3] (c) have intermediate radii (0.3 < r < 
0.7) between radii of LS and HS states. They feel high 



stress and the spin states of them are changing. 

In Fig. [4j we plot as triangles AE for LS domains at 
the side as a function of N^ s . Unlike the case of the 
corner, AE is not a simple decreasing function of Nl s 
for larger N^ s . It is worth noting that N^ s is no more 
equal to Nl s for 40 < A^ s , and takes a value for the 
configuration in the stationary state, which is given in a 
parenthesis in Fig. [4j LS domains which are bigger than 
a critical size are unstable at the side. 

If we set a round LS domain in the center of the HS 
phase, due to a huge distortion, the domain becomes very 
unstable and it collapses easily to reduce the number of 
LS molecules. 

We checked the qualitatively same tendency for the 
stability of LS domains when the system size (L) is larger. 
These observations lead to a major conclusion; LS do- 
mains cannot grow from the sides or the inner part of 
the system, which is different from the results of short- 
range interaction models. 

This conclusion suggests that different nucleation pro- 
cesses exist between short-range interaction models and 
the elastic model. We define Pi (corner), P2 (side), and 
P3 (inner part) as the nucleation rate from a corner, a 
side, and the inner part, respectively. Generally, they are 
functions of the size and temperature, and P\ > P2 > P3 
because the surface energy increases in this order. How- 
ever, if we take into account the possible location of nucle- 
ation, the probabilities to observe nucleation at a corner, 
a side, and the inner part of the system are Pi, L x P 2 , 
and L 2 x P 3 , respectively. Then, in short-range models, 
the relation Pi < L x P 2 < L 2 x P 3 will hold for larger 
L (linear dimension), and nucleation occurs in the inner 
part in large systems. Thus, so called multi-nucleation 
process takes place [29| . 

In the elastic model, however, P2 and P3 are essen- 
tially zero and Pi (corner) is the only probability for 
nucleation. Therefore, even if the system size is large, 
nucleation (clustering) always starts from corners. 

We next investigate the process in heating. Snapshots 
of transient states from the LS to HS phase in the heat- 
ing are given in Fig. [5] (a) — (d). Here, we also find local 
clusters of HS molecules around the corners, but in con- 
trast to the case of the process from the HS to LS phase 
(left branch of the hysteresis loop), a large homogeneous 
region appears as is observed in periodic boundary con- 
ditions. 

We consider the reason for the difference of the chang- 
ing pattern between the heating (LS to HS) and cooling 
(HS to LS) processes. In the SC system, the HS and 
LS states are not equivalent and we may expect different 
types of relaxations for the cooling and heating processes. 
In the cooling process at low temperatures, the energy 
stability is more important than the entropy gain and 
the nucleation from a corner is the most favorable. In 
the heating process at high temperatures, however, the 
entropy gain becomes more important, and the configu- 
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FIG. 5: (color online) Snapshots of configurations in the 
heating process from the LS to HS phase. Red (blue) circles 
denote HS (LS) molecules, (a) T = 1.0, t = 5452, (b) T = 1.0, 
t = 5522, (c) T = 1.0, t = 5672, (d) T = 1.0, t = 5748. 



ration may change uniformly, which can be seen in the 
inner part of the system. 

In summary, we have studied effects of the boundary 
condition in a SC model with effective long-range inter- 
actions. We found that domains always grow from cor- 
ners, which exhibits a striking contrast to the cases of 
short-range interaction models. In the heating process, 
an entropy-driven mechanism causes a smearing of clus- 
ters, and the configuration is close to that with the peri- 
odic boundary condition. 

The existence of macroscopic domains in SC com- 
pounds has been suggested in experimental studies of X- 
ray diffraction [30|, EE EH • The present study could give 
an insight into that suggestion. Dynamical properties 
of SC materials with OBC are important for studies of 
nano-scale systems, where the boundary plays a crucial 
role. 
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